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CNj Abstract. We present the first separate spectra of both components of the small-separation double QSO HE 0512-3329 

obtained with HST/STIS in the optical and near UV. The similarities especially of the emission line profiles and redshifts 
' strongly suggest that this system really consists of two lensed images of one and the same source. 

The emission line flux ratios are assumed to be unaffected by microlensing and are used to study the differential extinction 
effects caused by the lensing galaxy. Fits of empirical laws show that the extinction properties seem to be different on both lines 
of sight. With our new results, HE 0512-3329 becomes one of the few extragalactic systems which show the 2175 A absorption 
feature, although the detection is only marginal. 

We then correct the continuum flux ratio for extinction to obtain the differential microlensing signal. Since this may still be 
^/-^ ' significantly affected by variability and time-delay effects, no detailled analysis of the microlensing is possible at the moment. 

This is the first time that differential extinction and microlensing could be separated unambiguously. We show that, at least in 
**"~^ HE 0512-3329, both effects contribute significantly to the spectral differences and one cannot be analysed without taking into 

account the other. For lens modelling purposes, the flux ratios can only be used after correcting for both effects. 
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. . , 1 . Introduction because of the wealth of absorption features and the close im- 

■ age separation. 
! > ■ HE 0512-3329 was discovered as a probable lensed quasar _ _ . . .... 

^r 1 . . , „ , . , , _ _ , The propagation effects of microlensing and extinction in 

r > , in the course of a snapshot survey with the Space Telescope , , . , , ,• , , , r 

5-H T „ . ,„_„. T . , ... , ^r,^ ., the lensing galaxy can be studied by using the fact that we see 

rrj Imaging Spectrograph (STIS). It is a doubly imaged QSO with . , , * . , . 

, . . ° , . . . , two images of the same source. In singly imaged sources, it is 

source redshitt of z — 1.58 and an image separation of ?. . . . . , ... . . 

not possible to uniquely identify microlensing or extinction et- 



m 

CO 
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0'.'644. The lensing galaxy has not been detected yet, but strong 
metal absorption lines with a redshift of z — 0.93 identified in 
the integrated spectrum provide good evidence that a damped 

Lye (DLA) system intervenes at this redshift. Th is system is , . ,. „, ..... . , 

... , j„ ri ual images as diagnostic. These differences are independent of 

very probably associated with the lensing galaxy (Gregg et al. , , _ , , . , 

" " the true source spectrum and flux and therefore allow the ex- 

amination of differential propagation effects. 

An additional complication arises from possible variability 



fects because the true source fluxes and spectra are not a priori 
known. In the case of multiply imaged sources, this problem 
can be circumvented by using only differences of the individ- 



2000). 

The absorption lines provide the possibility to study the 

interstellar matter (ISM) of the lensing galaxy along the two 

,. c ■ * j u e 1 1 . j i j of the source. Because of the time-delay, observations of both 

lines of sight, separated by 5.1kpc at the proposed lens red- ■" 

, l . • . ii .1 u . t components taken at the same time are really snapshots of the 

shift, which allows us to learn more about the nature of the r ; 

tit * u u u • i j ex. j j source at different epochs. To correct for this effect (which is 

DLA absorber by comparing column densities of hydrogen and r v 

, , ii i • £i i u l- c ■ thought to be negligible for the emission lines because of the 

metals as well as kinematic profiles along both lines of sight. b b b 

Ilr rvi-n TT-,n • • n n . j£ u- " longer variability time scale), (spectro)photometnc monitoring 

HE 0512-3329 is especially well suited for such investigations ° , . , , 

is necessary to determine the time-delay and obtain measure- 



Send offprint requests to: O. Wucknitz ments at the same source e P och ' 

1 We assume a low-density flat universe with Q. = 0.3, A = 0.7 and It has been known from earlier photometric observations 
H = 70kms~' Mpc~' throughout this paper. that the flux ratio A/B in HE 0512-3329 shows a strong de- 
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pendence on wavelength. In the R and / bands, A is brighter 
than B by about 0.45 m ag while the two are almost equal 
in B iGre^.g_ej "aL 2000). A natural explanation for this effect 
is differential reddening caused by different extinction effects 
in the two lines of sight. This scenario is supported by the 
strong metal lin es already visible in the integrated spectrum of 
HE 0512-3329 l lGregg et al. 2000) which make the existence 
of signifi cant amounts o f dust very likely. Under this assump- 
tion. lOregg et alJ J2000l) estimated a relative visual extinction 
of Ay - 0.34 mag by fitting galactic extinction curves with 
Rv = E(B-V)/Av = 3.1 on both lines of sight to the photomet- 
ric data. We will see later that this approach is not appropriate 
in the case of HE 0512-3329. 

The other possible explanation is microlensing by stars and 
other compact objects in the lensing galaxy. Since the lines of 
sight cross different star fields in the lens, the microlensing ef- 
fect is independent for A and B and can be observed in the dif- 
ference. The lensing effect itself is an achromatic process, but 
chromatic effects are nevertheless expected if different parts of 
the source, which are influenced by different microlensing am- 
plifications, have different colours. Since the microlensing am- 
plification pattern has a high contrast on very small scales, the 
effectively observed total amplification of the integrated source 
depends significantly on the size of the source. The effect is 
maximal for point sources but is smoothed out and thus reduced 
for larger sources. 

Standard models of QSOs consist of different emitting re- 
gions of very different sizes (see e.g. lKrolildll999t) . Smallest 
is the continuum region, which itself scales with wavelength. 
Broad line regions are expected to be much larger than the 
continuum region but smaller than the narrow line regions. 
The whole range of scale lengths covers several orders of 
magnitude. This means that the microlensing signal should be 
strongest for the continuum, only very weak at most for the 
broad emission lines and absent for narrow lines. 

This size-dependent behaviour of microlensing gives a pos- 
sible clue to separate the effect from differential extinction and 
study the two processes separately. Microlensing is smeared 
out in the emission lines whose fluxes are thus only influ- 
enced by the differential extinction. The continuum, on the 
other hand, is also affected by microlensing. The extinction it- 
self is a smooth function of wavelength but does not depend 
on the source size. We can therefore use the emission lines to 
study the extinction directly without worrying about microlens- 
ing and correct the continuum for extinction to study the pure 
differential microlensing in the remaining differences. 

The spectra we use to follow this approach were obtained 
with STIS on board the HST in the optical and near UV cover- 
ing the QSO emission lines of Ly B, Ly a, Si iv, C iv and C m] 
with the goal to study the differential extinction by comparing 
the continua of the A and B images. 

The results of these observations are presented in this paper. 
We show the first separate spectra of the A and B images and 
use the similarity of these to confirm the nature of the object as 
a gravitational lens system. We use the combination of contin- 
uum and emission lines to disentangle the effects of extinction 
and microlensing as described above. We see that microlensing 



actually does influence the flux ratio significantly and that the 
continuum cannot be used to study the extinction directly. 

To analyse the differential extinction we will fit empirical 
extinction laws to our emission line data and discuss the results. 
The microlensing data will be used to estimate an upper limit of 
the source size using the assumption that variability and time- 
delay effects can be neglected. 

This is the first time that such a detailed attempt to separate 
the effects of differential extinction and microlensing has been 
successfully performed for any lens system. 

2. Observations and initial data reduction 

To cover a sufficiently large spectral range, we used the CCD 
for the optical and the near UV MAMA detector for the UV 
part. The slit was positioned along the line connecting the two 
QSO components. 

The CCD spectra were taken on 13 August 2001 with the 
G430L grating and 52x0.5 apert ure, coverin g a spectral range 
from 2900 to 5700 A JLeitherer etalJl200lh . With a spectral 
element size of 2.73 A, the resolution is about 3.8 A. Three ex- 
posures with total integration time of 2266 sec and a standard 
dithering pattern were used. 

Data reduction started with a standard 2-dimensional re- 
duction including flat-fielding, flux and wa velength calibration 
with the CALSTIS software Version 2.10 ( iBrown et al] Eo02) 
in IRAF V2.12 with STSDAS Version 2.3. Cosmic ray rejec- 
tion could not be performed in a standard way because the indi- 
vidual exposures were not done in CR-SPLIT mode. We there- 
fore used our own software, developed specifically for this task, 
to apply the shifts, reject cosmic ray events and combine the 
images. 

We then used the CALSTIS task X1D to extract the spectra 
of both QSO components separately. A smaller than standard 
extraction width was used to minimize the fraction of measure- 
ments being rejected due to flagged bad pixels or cosmic rays. 

The three resulting extracted spectra were finally combined 
and cleaned of still remaining cosmic rays and otherwise de- 
viant pixels with our own software. 

The MAMA spectra were taken with the G230L grating 
and 52x0.2 aperture, covering a range from 1570 to 3180A 
with an element size of 1.58 A and a resolution of 3.3 A. The 
six exposures of a total of 16588 sec were taken on 15 August 
2001, only two days after the CCD spectra. Since the MAMA 
detector is not sensitive to cosmics, a standard CALSTIS re- 
duction could be used to calibrate the 2-dimensional images. 
To extract the two QSO components separately, X1D was used 
with a slightly smaller than standard extraction width and with 
background regions adapted to the situation of two closely sep- 
arated but well resolved spectra. The six MAMA spectra were 
then combined with the same algorithm as the CCD spectra. 

For the analysis, we did not use the short wavelength end 
of the combined CCD spectrum (< 3020 A) and the long wave- 
length end of the MAMA spectrum (> 3 120 A). In these re- 
gions, the S IN is not optimal and calibration uncertainties are 
expected to be worst. In the remaining overlapping region, the 
two spectra agree very well and a simple averaging combina- 
tion was used for the final spectrum. 
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Fig. 1. The combined spectra of HE 0512-3329, smoothed with a Gaussian filter of FWHM=3.5 A. Component A is at top, B at 
bottom. The curve near fx — shows the statistical errors. The middle panel marks emission and absorption features at z — 1 .58 
(QSO emission), z = 0.933 (DLA system), z — 1.137 (Lyman limit system) and z — (galactic absorption). 
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The typical error per spectral element of the final com- 
bined spectrum is 1-2 • 10~ 17 erg s _1 ctrT 2 A -1 , except near 
A ~ 3 100 A and for A < 1700 A. This is equivalent to a S/N of 
30 per element in the optical and about 15 in the UV region. 



3. The spectra 

3.1. First view 

Figure ^ shows the combined spectra from all CCD and 
MAMA exposures. These comprise the first separate spectra 
of the two QSO components in HE 0512-3329. The general 
features of A and B look very similar. The relative strength of 
the emission lines, their profiles (see also Fig.|4]and discussion 
below) and redshifts agree very well which strongly confirms 
the interpretation as two lensed images of one source. The con- 
tinua, on the other hand, show similarities as well as significant 
differences which will be discussed in detail in the later parts 
of this paper. 

The wealth of absorption features seen in both images can 
provide important information about column densities of both 
hydrogen and metals and will be discussed in detail in another 
publication. The spectra confirm the existence of a DLA sys- 
tem at a red s hift of z = 0.933 which was already proposed by 
ICreggetalJ lfcoOO) to explain the strong low-ionization metal 
lines in the (then still combined A+B) spectrum. Our estimate 
of the r edshift differs only slightly from the former measure- 
ment by |Gregietal] l|200(l of z = 0.9313 + 0.0005. Figure □ 
identifies the anticipated DLA line as well as a high number of 
metal absorption features at the same redshift in both compo- 
nents. The break at about 2000 A, which is especially promi- 
nent in B due to the higher continuum there, does not corre- 
spond to this already known system at z — 0.93 but is caused 
by a second intervening Lyman limit system at a redshift of 
about z = 1.137. This system is much weaker than the lower 
redshift one but nevertheless shows some metal lines which are 
also identified in Fig.[Q Weak Mgn absorption at a very sim- 
ilar re dshift (z — 1.1346) was already noticed bv lGregg et alJ 
(2000). Finally we detect galactic Mg i+n absorption at a red- 
shift consistent with zero. 

A quantitative analysis of the absorption features, com- 
bined with high resolution UVES spectra, will be published 
in an upcoming paper (Lopez et al., in prep.). Especially im- 
portant is the analysis of the DLA lines in our spectra which 
provide the H i column density needed to determine metalici- 
ties in both lines of sight. 

3.2. Differences in the global shape of A and B 

Although the two spectra have many features in common, there 
are some highly significant differences. Most prominent is the 
very different behaviour of the continuum flux on large wave- 
length scales. On the red side of A « 4600 A, the A component 
is brighter than B. For smaller wavelengths, especially close to 
the limit near 2000 A, B becomes much brighter than A . This 
effect has already been described by Greg g et all d20 00) using 
photometric data in B, V, R and /. While A is brighter by about 



0.45 mag in / and R, the difference starts to decrease in V and 
is compatible with zero in the B band. 

Given the similarities of the spectra, the differences have to 
be interpreted in a lensing context. As discussed in the intro- 
duction, two possible explanations will be discussed here. One 
is differential reddening by different amounts of dust along both 
lines of sight in the lensing galaxy, the other is microlensing 
which by different amplification patterns for the two images in 
combination with a colour gradient of the source can also result 
in chromatic effects. 

To separate these effects, we use the assumption that the re- 
gions causing the emission lines are much larger than the con- 
tinuum region so that the microlensing signal is smeared out 
in the emission lines. We can thus use the emission lines to 
study the pure extinction signal. In the next step we correct the 
continuum (which is affected by extinction and microlensing) 
for extinction to study the microlensing effect. To accomplish 
these tasks, we have to fit the continuum and determine the flux 
ratios of the continuum-subtracted emission lines. 

3.3. Global continuum fit 

We first flagged the regions of the five strong emission lines of 
both spectra and kept only the continuum which is still covered 
with many absorption lines. We tried two approaches to fit the 
global continuum. One is a global fit of a polynomial of degree 
4 to the logarithmic flux. The other is a smoothing algorithm 
which works by fitting a quadratic function to the spectrum for 
each wavelength, using a Gaussian weighting with FWHM of 
300 A to effectively include only nearby data points. The cen- 
tral value of this local fit is taken as the result for the central 
wavelength and the procedure is repeated for each wavelength. 
This algorithm is more appropriate than a simple Gaussian con- 
volution filter because it can follow local slopes and curvatures 
more accurately. In particular it avoids any edge effects. In both 
algorithms we allowed for a certain fraction of pixels to devi- 
ate in a negative direction from the model to avoid absorption 
lines affecting the fit. These pixels were then discarded in the 
X 2 summation. We started with this fraction set to zero and then 
increased it slowly until a level of 30 % was reached. This value 
was chosen to obtain a final normalized x 1 close to unity. It is 
equivalent to a final effective (one-sided) clipping limit of about 
1.5-1.7(7-. 

The global continuum fits are shown in Fig. |2] together 
with the observed spectrum. In the comparison, the locally 
smoothed continuum follows the spectrum more closely while 
the global polynomial fit is slightly smoother on large scales. 
Both curves look like good continuum fits to the eye, and the 
small differences are not significant for the further analysis 
which was performed with both fits. We have to keep in mind, 
however, that the continuum short-wards of the Ly B emission 
line might still be significantly contaminated by the Lyman for- 
est. This should not be a problem in the continuum subtraction 
of the emission line, because the broad emission lines are af- 
fected by the Lyman forest as well. The shape of the continuum 
near the short wavelength end itself should nevertheless be in- 
terpreted with caution. Figure|3]shows the continuum of A and 
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Fig. 2. Observed spectrum of A (top) and B (bottom) together with the two different continuum fits. The global polynomial fit 
is shown solid, the locally smoothed estimate is dashed. The shaded regions mark the emission lines and the region affected by 
the Lyman limit on the far left. Parts of the spectra not included in the fit (in the emission line regions or rejected by the global 
polynomial fitting algorithm because of low flux) are shown dotted, the rest of the spectra is solid. 



B (top) as well as the ratio of the two (bottom) as a function of 
A. 

3.4. Emission line fluxes 

The main problem in calculating the emission line fluxes is the 
continuum subtraction. We used three different estimates for 
the continuum to get a feeling for the accuracy. The first two 
are the global polynomial fit and the smoothed continuum as 
discussed before. The third method uses local linear trends with 
fitting regions on both sides and directly next to the emission 
lines. After subtracting the continuum, different methods were 
used to measure the emission line flux. The simplest one uses 
weighted averaging over two different apertures, one narrow 
(27-59 A) and one relatively wide (74-261 A). The widths of 
these windows were adapted 'by eye' to the individual emis- 
sion lines. The narrow windows just enclose the core of the 
lines while the wide ones are chosen in a way to include the 
complete emission lines. We also performed Gaussian fits and 
a local quadratical smoothing (as explained before, Gaussian 
FWHM=100 A) of the flux, also with wide and narrow fitting 



regions. The values at the nominal line centre of these fits were 
then used as the maximal line flux. 

A byproduct of the Gaussian fits is the QSO redshift. For 
the five emission lines, we obtain a mean value of z = 1.582 
with a rms scatter of 0.008 which is dominated by the differ- 
ent line shapes. The statistical error is much smaller, of the or- 
der 0.001. This com pares moderately well with the result from 
Greg g et ai1 d2000) for the combined spectrum of z = 1 -565 for 
C iv and C m] . There is no noticeable difference in the results 
for A and B. 

Figure [5] (bottom) shows the resulting emission line flux 
ratios compared with the continuum ratio. The plotted values 
are the mean results from our different estimates. As error bars 
we used the formal statistical error added in quadrature to the 
internal scatter of all methods to include possible systematics. 
There seems to be almost no variation for the lines C iv, Si iv, 
Ly a and Ly ft. All of these are compatible with their mean of 
log(/A//B) = 0.04. Only for Cm] we notice a deviation from 
the other valu es. The photometric BVRI measurements in that 
plot are from iGregg et alJ (|2000), H is from a deep 2400 sec 
VLT/ISAAC exposure (details will be published elsewhere). 
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Fig. 3. Top: Continuum flux of A and B. Bottom: Flux ratio A/B (curve). The error crosses above the lines are for the continuum 
subtracted emission lines, where the horizontal bar shows the FWHM of the lines. Photometric measurements (combined contin- 
uum and emission lines) are shown at the right (see text). The photometric error bars are formal statistical uncertainties. R, I and 
H are outside of the wavelength range and are shown as arrows. The same line type coding as in Fig.|2]is used for the continuum 
fits in both panels. The regions of emission lines were not used for the fits and are shown dotted. 



Figure |4] shows the emission line profiles after subtraction 
of the local linear continuum fit and rescaling with the esti- 
mated emission line flux ratio. Apart from subtle differences 
in absorption features, the emission line profiles of A and B 
are remarkably similar. The remaining differences seem to be 
largest in Si iv but this is a result of the weakness of this line 
which also leads to the relatively large photometric uncertainty 
shown in Fig. [3] 

The differences of the absorption left of Ly ft are a result of 
the very different relative strength of the absorption systems in 
A and B. Altogether the profiles are highly consistent, leaving 
very little doubt about the lensed nature of the two images. 

The emission line ratios determined from the wide (com- 
plete line) and narrow (line core) windows are also compatible 
with each other. This means that the microlensing signal is the 
same for the broad and narrow line regions. Since the two are 
expected to be very different in size, the effect would be dif- 
ferent if there is any effect at all. This means that microlensing 
seems to be absent already for the broad lines. Using emis- 
sion lines therefore seems a reasonable way to measure extinc- 
tion curves which are unaffected by microlensing, at least in 
HE 0512-3329. 

4. Differential extinction 

To analyse the differential extinction, we fitted different em- 
pirical extinction laws to the emission line data and the photo- 
metric data in / and H, assuming that microlensing is affect- 
ing these bands only weakly as indicated from the similarity 
of the extrapolated continuum and the photometric data (Fig.|3| 
bottom). To take into account the residual microlensing effects 
and flux ratio changes by variability within the time-delay, we 



used an error of <x = 0.05 mag for the photometric data points 
instead of the formal error or 0.02 mag. The results are not very 
sensitive to this exact value. 

The first e xtincti on law is the galactic curve from 
ICardelli et all Jl989l) which has a 2175 A feature 
whose strength is correlated to the far-UV behaviour 
(CCM). Alternative l y we used the parametrisation from 
Fitz patrick & Massal {l990) once with the Drude term de- 
scribing the 2175 A bump (FM bump) and once without it 
(FM nobump). The strength of the bump is a free parameter 
independent of the global shape of the extinction law, but we 
fixed its position and width {xq = 4.6 u rn" 1 , y = Q.95 urn" 1 ). 
We n eglected the far-UV term C4F(x) (Fitz patrick & M assa 
1990, eq. 2) because it becomes important only at shorter 
wavelengths and is thus only very weakly constrained. The 
inclusion would actually increase the normalized x 2 of the 
fits. The strength of the 2175 A feature does not depend on 
the far-UV term. Although the FM models are valid only 
for /I -1 > 3.3 |xm~ , we included the H and / data in the fit 
using an extrapolation. The FM plots should therefore only be 
seen as illustrations and not as physical models. Without H 
and / these curves would rise on the left to match C m] more 
closely, without changing the 2175 A bump very much. As a 
second extragalacti c model we used the functional form from 
ICalzetti et ail fcOOOh (CAB) without a 2175 A feature. All 
extinction law fits are shown in Fig. [3] 

We see that all of these models fit the data relatively well; 
the normalized \ 2 values are all close to unity. To follow the 
slightly deviating value from the C iv line more closely, how- 
ever, the 2175 A bump has to be included which also shows in 
improved residuals for the CCM and FM bump models. Since 
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Fig. 4. Profiles of the five emission lines after local continuum subtraction. A is plotted solid, B is dashed. The flux of B was 
scaled by the determined emission line flux ratio to obtain the same level for both components. The shaded regions were used for 
the linear continuum fit with rejection of absorption features. The spectra are smoothed with a Gaussian kernel of FWHM=5 A. 
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Fig. 5. Different extinction laws fitted to the data from emission lines and red/infrared broad band photometry. The photometric 
error bars were extended to ±0.05 mag (see text). Models with 2175 A bump: solid (FM bump) and dotted (CCM); models 
without: dashed (CAB) and dot-dashed (FM nobump). 



the C iv line lies in a region of the spectrum were the contin- 
uum is very well defined, there is no reason to trust this data 
point less than the others. Our result can therefore be seen as 
a marginal detection of the 2175 A feature in the spectrum of 
HE 05 12-3329. This interpretation is supported by the fact that 
the free fit of the absorption feature in the FM bump model 
favours a strength almost exactly equal to the galactic CCM 
model. 



For the physical interpretation we only consider the galactic 
CCM law 



MA) 
A(V) 



= a(A~ l ) + 



Rv 



(1) 



with Ry = A(V)/E(B - V). The functions a and b are de- 
fined piecewise over a region 0.3 uirT 1 < /T 1 < lOixtrT 1 
dCardelli et aljfl989i) . For the differential extinction law for 
A - B we obtain 



A\X) - A B (A) = p a a(A- 1 ) + p h b(A~ l ) 



(2) 
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Table 1. Parameters of CCM model fits in magnitudes, po is 
the 'true' brightness difference, corrected for extinction but in- 
cluding the different lens amplifications. p a is the differential 
extinction in V, pi, is the differential colour excess. The effec- 
tive R v is calculated as p a j pt,. It is the value of R v if the shape 
of the extinction laws in A and B would be equal or if extinc- 
tion is affecting only one component. The first two lines show 
a free fit including 1 <x error bars, R e y was fixed for the other 
lines. 



Po 


Pa 


Pb 


Rf 


x 2 /y 


-0.493 


-0.074 


+0.037 


-2.00 


3.2/4 


±0.11 


±0.11 


±0.015 




(free fit) 


-0.672 


+0.131 


+0.059 


+2.20 


6.1/5 


-0.733 


+0.206 


+0.067 


+3.10 


8.8/5 


-0.870 


+0.447 


+0.077 


+5.80 


34.3/5 




Fig. 6. Residuals for fixed parameters p a and pi, of the CCM 
model. The ellipses show differences of x 1 from its minimal 
value of 1, 4 and 9 corresponding to confidence limits of 1, 2 
and 3 cr for one parameter. The solid lines intersecting at (0, 0) 
show the effective Ry = p a lpb- We see that the data are com- 
patible with one value of Ry only if this value is very low. The 
best fit has p a /pb = -2.00. The dotted parallel lines show the 
fitted po which gives the true amplification ratio of the images. 



with the two independent parameters 
Pa = A A (V)-A B (V) 



Pb = 



A A (V) A B (V) 



(3) 
(4) 



The first parameter p a is the differential visual extinction, the 
second one is the differential colour excess ph = E A (B — V) — 
E B (B - V). The difference of measured magnitudes is equal to 
the difference in intrinsic brightness (caused by different ampli- 
fication and variability in combination with the time-delay) po 
plus the difference in extinction. We therefore have three free 
parameters. 

The parameters of the best fit are shown in Table ^ We 
see that the visual extinction A(V) has to be higher in B by 



0.07 mag but the reddening E(B-V) is higher in A by 0.04 mag. 
Recall that A(V) and E(B - V) are not measured directly in V 
and B but are merely used as parameters to describe and fit the 
UV extinction law. Figure |6] shows residuals of fits with dif- 
ferent fixed values of p a and pt,. The data are only marginally 
compatible with only one value of Ry (same value in both im- 
ages or extinction in only one of the images) in the realistic 
region of 2.2 < Ry < 5.8. The best effective Ry - p„lpb 
is -2.00 which can only be explained with a combination of 
extinction both in A and B with different Ry. To explain the 
measurements with the smallest possible visual extinction, Ry 
must be at its lower and R B at its upper limit. For Ry = 2.2 and 
R B = 5.8 we get A A (V) = 0.177 mag and A B (V) = 0.251 mag. 
For more similar extinction laws, the total extinction has to be 
higher, e.g. A A (V) = 0.668 mag and A B (V) = 0.742 mag for 
R A = 3 and R B = 4. Such values are not atypical for galactic 
extinction. Even with more similar Ry the total extinction does 
not become unrealistically high. 

The highly significant differences between continuum and 
emission line ratios (Fig.|3|l show the importance of microlens- 
ing when studying differential extinction curves. Using photo- 
metric data points without correcting for microlensing like in 
iFalco et alJ dl999h would in the case of HE 05 12-3329 lead to 
highly incorrect results. 

An important result for lens models is the true amplification 
ratio which corresponds to po = A-B = -0.49 + 0.1 1 mag (1 cr 
errorbar). The amplification ratio is thus |jUa///b| = 1.57 + 0.17. 
Note how different this value is from the continuum or emission 
line ratio in the UV (Fig. The combination of extinction 
and microlensing changes the continuum flux ratio at the short 
wavelength end by almost 2 magnitudes. One should therefore 
be very careful when using optical flux ratios as constraints for 
the amplification ratio in lens modelling if no correction for 
microlensing and differential extinction can be applied. The re- 
sult for \pa/p~b\ depends on the emission lines only. Variations 
of the continuum during the time-delay (see estimate below) 
would not change this ratio. Variations of the emission lines 
are still a possible source of error here but are expected to be 
much smaller because of the much larger size of the emitting 
regions. 



5. The lensing galaxy 

The available observational data are quite insufficient to con- 
strain realistic lens models to a satisfying degree. Since not 
even an estimate of the galaxy position is available, the best 
we can do at the moment is to use the simplest model, a singu- 
lar isothermal sphere without external perturbations, and esti- 
mate the mass with this approach. In contrast to the time-delay, 
the mass depends only weakly on the exact galaxy position or 
asymmetries of the mass distribution. 

With the source redshift z s = 1.58 and a proposed lens 
redshift Zd = 0.933, this leads to an expected radial veloc- 
ity dispersion of cr v = 192kms~ 1 and an enclosed mass 
within the Einstein radius of M E = 6.8 • 1O 1O M0. Using the 
Tully -Fisher relation for late-type galaxies at intermediate red- 
shift (Zieeler et alJl2002l) . we obtain an absolute luminosity 
of M B = -22 mag using an equivalent circular velocity of 
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v c = 330km s . This corresponds to an apparent / magnitude 
of 21.5 mag, about 4.5 mag weaker than each QSO component. 
Direct images taken with STIS in the same campaign do not 
show any obvious sign of the core of the lensing galaxy after 
applying a standard subtraction of the two QSO images. We 
are currently working on new methods to optimize this PSF 
subtraction process. 

The maximal time-delay would be reached if the lens is 
located very close to one of the images, leading to a limiting 
value of 170 days. For this scenario a very high ellipticity of 
the mass distribution would be required. More realistic values 
are ~ 1-2 months. A spherically symmetric lens with a true am- 
plification ratio of 1 .57 (the result of the previous section and 
assuming negligible variability) would lead to a time-delay of 
Af = 38 days with A leading. Even a moderate ellipticity or 
external shear could change this number (even its sign) consid- 
erably. 



6. Microlensing 

With the measurement of the emission line flux ratios, we were 
able to obtain estimates of the differential extinction. Since the 
continuum is affected by extinction and microlensing, we have 
to remove the extinction effect to study the microlensing signal 
alone. Figure shows the continuum flux ratio corrected for 
differential extinction as determined from the emission lines. 
Given the high uncertainties in estimating the extinction law, 
we do not show the estimated global microlensing function but 
only the values for the wavelengths of the emission lines which 
have been measured directly and which are thus not model de- 
pendent. 

These data points should be flat zero in the absence of mi- 
crolensing. We see that in reality either B is amplified and/or 
A is de-amplified with a strength increasing for smaller wave- 
lengths. 

An additional complication is the possibility of variabil- 
ity. Because of the time-delay, the two components are ob- 
served in two different epochs in the source plane. If the source 
brightness or spectrum changes significantly over this time- 
delay, the measured ratio would not reflect the true ratio any- 
more. Assuming that this effect is negligible in our observa- 
tions of HE 0512-3329, we can use the magnitude of the 
differential microlensing to estimate an upper limit for the 
sour ce size in terms of the Einstein radius using the recipe 
of lRefsdal&^t abell (1991). If the observed differential mi- 
crolensing signal of 5m = 0.5 mag is typical for the mi- 
crolensing variance, a conservatively estimated upper limit of 
4 Einstein radii can be derived. The Einstein radius in our 
case is 8 e — 1 .2 u.arcsec ■ ^JM/Mq for microlenses of mass 
M, corresponding to 0.01 pc (2200 AU) in the source plane for 
M = M Q . 

We also tried to fit the parameters of more specific physical 
scenarios like isolated caustic crossings of a standard accretion 
disk model to the data, only to learn that (although this should 
be a powerful method in principle) the available data are not 
sufficient to obtain any conclusive results. 
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Fig. 7. Differential microlensing amplification (A/B) as a func- 
tion of wavelength (error bars). The continuum (affected by mi- 
crolensing and extinction) was corrected for extinction using 
the emission line fluxes which are expected to be unaffected by 
microlensing. 

7. Summary 

We presented the first separate spectra of both components of 
the doubly lensed QSO HE 0512-3329. The similarity of the 
optical and near UV spectra strongly suggest the lens nature of 
this object. 

The known optical colour differences of the images, which 
(together with the metal absorption lines) were the main moti- 
vation to start this project, extend to the near UV and become 
even stronger in this range. While A is brighter than B by about 
0.5 mag in R and / band, this relation is more than inverted 
in the continuum near 2000 A where B is brighter by about 
1.3 mag. 

We used the continuum subtracted emission line flux ra- 
tios to study the differential extinction in both images. The 
line emitting regions of the QSO are confirmed to be so large 
that any residual microlensing effects on their fluxes is in- 
significant so that the extinction can be studied unambiguously. 
We fitted several empirical extinction laws from the literature 
and f ound that curves with a signific ant 2175 A bump (C CM: 
ICardelli et all l989: FM bump: iFitzpatrick & Mas"salll990l) re- 
produce the data better than typical extragalactic curves with- 
out this feature (FM nobump: Fitz patrick & Massall990t CAB: 
Calz etti et al]l2000n . Although the differences are not highly 
significant, this provides at least evidence that the 2175 A fea- 
ture can also be present in high redshift galaxies (compare 
iMotta et all2002l 

From the difference between the emission line and contin- 
uum ratios we learned that microlensing contributes signifi- 
cantly to the variations of flux ratio with wavelength. On the 
short wavelength end, the microlensing signal is actually much 
stronger than the differential reddening. While the emission 
lines have almost the same flux in A and B there, the continuum 
is brighter in B by up to 1.3 mag (see Fig. [3] bottom). For our 
study it was therefore absolutely essential to use the emission 
line fluxes for the extinction analys is. Analysing the c ontinuum 
or broad band photometry like in iFalco et all Jl999h without 
correcting for microlensing would have led to highly incorrect 
results in the case of HE 0512-3329. 

The analysis of galactic extinction law fits showed that the 
data are only marginally consistent with the same extinction 
law in A and B with only different strength. The best fit has a 
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higher total visual extinction A(V) in B but a higher reddening 
E(B - V) in A. 

Using the image separation, we were able to estimate the 
mass within the Einstein radius to be Me = 6.8 • 10 1() Mo- 
The expected radial velocity dispersion of an isothermal model 
is <j v = 192kms~ I , equivalent to a circular velocity of v c - 
330 km s -1 . The time-delay is expected to be of the order 1-2 
months. 

In the last part we used the differential extinction measure- 
ment from the emission lines to correct the continuum flux ratio 
at five wavelengths for this effect. The remaining ratio is a di- 
rect measurement of the microlensing signal alone. This is the 
first time that such a study is carried out to separate the two 
effects and discuss them individually. 

The microlensing data are not sufficient to determine the 
parameters of physically motivated source models accurately 
but could be used to estimate an upper limit for the effective 
source size using the assumption that variability over a scale of 
the time-delay does not change the observed flux ratio signifi- 
cantly. To test this assumption and to determine the time-delay, 
photometric monitoring of this system should be performed in 
the future. The information gained by our work (spectral mi- 
crolensing measurements at one epoch) is complementary to 
the more common light-curve analysis techniqu e (photometric 
measurements at many epochs, see e.g. Shalvapin etal]l2002l 
for recent results). A combination of both methods by using 
spectrophotometric monitoring data has the potential to lead to 
more valuable constraints than both methods alone by break- 
ing degeneracies inherent in the individual approaches. It will 
then be possible to determine brightness profiles in all of the 
measured wavelengths which can be combined to a relatively 
model-independent measurement of the source's temperature 
profile, allowing to test different physical QSO models. 
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